Rubidium and lead abundances in giant stars of the globular 

clusters M 13 and NGC 6752^ 



David Yong 

Department of Physics & Astronomy, University of North Carolina, Chapel Hill, NC 

27599-3255 

yongOphysics . unc . edu 

Wako Aoki 

National Astronomical Observatory, Mitaka, 181-8588 Tokyo, Japan 

aoki . wakoOnao .ac.jp 
David L. Lambert 
Department of Astronomy, University of Texas, Austin, TX 18112 
dllQastro . as . utexas . edu 
Diane B. Paulson 

NASA's Goddard Space Flight Center, Code 693.0, Creenhelt MD 20111 
diane . b . paulson@gsf c . nasa . gov 

ABSTRACT 

We present measurements of the neutron-capture elements Rb and Pb in 
five giant stars of the globular cluster NGC 6752 and Pb measurements in four 
giants of the globular cluster M 13. The abundances were derived by comparing 
synthetic spectra with high resolution, high signal-to-noise ratio spectra obtained 
using HDS on the Subaru telescope and MIKE on the Magellan telescope. The 
program stars span the range of the 0-Al abundance variation. In NGC 6752, 
the mean abundances are [Rb/Fe] = -0.17 ± 0.06 {a = 0.14), [Rb/Zr] = -0.12 
± 0.06 (a = 0.13), and [Pb/Fe] = -0.17 ± 0.04 (a = 0.08). In M 13 the mean 
abundance is [Pb/Fe] = —0.28 ± 0.03 (cr = 0.06). Within the measurement 
uncertainties, we find no evidence for a star-to-star variation for either Rb or 
Pb within these clusters. None of the abundance ratios [Rb/Fe], [Rb/Zr], or 
[Pb/Fe] are correlated with the Al abundance. NGC 6752 may have slightly 
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lower abundances of [Rb/Fe] and [Rb/Zr] compared to the small sample of field 
stars at the same metallicity. For M 13 and NGC 6752 the Pb abundances are in 
accord with predictions from a Galactic chemical evolution model. If metal-poor 
intermediate-mass asymptotic giant branch stars did produce the globular cluster 
abundance anomalies, then such stars do not synthesize significant quantities of 
Rb or Pb. Alternatively, if such stars do synthesize large amounts of Rb or 
Pb, then they are not responsible for the abundance anomalies seen in globular 
clusters. 

Subject headings: stars: abundances - Galaxy: abundances - globular clusters: 
individual (M 13, NGC 6752) 



1. Introduction 

Globular clusters are ideal laboratories for testing the predictions of stellar evolution the- 
ory (e.g., Renzini & Fusi Pecci 1988) since the individual stars arc believed to be monometal- 
lic, coeval, and at the same distance. In a given globular cluster (excluding uj Cen), spec- 
troscopic observations of individual stars have confirmed that members have uniform com- 
positions, at least for the Fe-peak elements (e.g., see review by Gratton ct al. 2004 and 
references therein). However, it has been known for many years now that globular clusters 
exhibit star-to-star abundance variations for the hght elements C, N, O, Na, Mg, and Al 
(e.g., see review by Kraft 1994). Specifically, the abundances of C and O are low when N 
is high and anticorrelations are found between O and Na as well as Mg and Al. Recently, 
variations in the abundance of fluorine have been discovered in giants in M 4 where the 
amphtude of the dispersion exceeds that of O (Smith et al. 2005). 

It is generally assumed that the light element variations arise from proton-capture re- 
actions (CNO-cycle, Ne-Na chain, and Mg-Al chain), though the specific nucleosynthetic 
site(s) remain elusive. One possibility for the origin of the star-to-star abundance variations 
is deep-mixing and internal nucleosynthesis within the observed stars. Evidence for this 
"evolutionary scenario" include C and N abundances (Suntzeff & Smith 1991) that vary 
with location on the red giant branch (RGB). Extensive mixing down to very hot layers is 
necessary to change the surface composition of Na, Mg, and Al. Such mixing is not predicted 



^ Based in part on data collected at the Subaru Telescope, which is operated by the National Astronom- 
ical Observatory of Japan and on observations made with the Magellan Clay Telescope at Las Campanas 
Observatory. 
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by standard models and the proposed mechanisms include meridional circulation (Sweigart 
& Mengel 1979), turbulent diffusion (Charbonnel 1995), and hydrogen-burning shell flashes 
(Fujimoto et al. 1999; Aikawa et al. 2001, 2004). An alternative possibility is pollution from 
intermediate-mass asymptotic giant branch stars (IM-AGBs), first suggested by Cottrell & 
Da Costa (1981) to explain the Na and Al enhancements observed in CN strong stars in 
NGC 6752. In IM-AGBs, hydrogen-burning at the base of the convective envelope, so-called 
hot bottom burning (HBB), can produce the observed C to Al abundance patterns. Either 
the ejecta from IM-AGBs pollute the proto-cluster gas from which the present cluster mem- 
bers form or the ejecta are accreted by present cluster members. The strongest evidence 
for this "primordial scenario" has come from observations of main-sequence stars in which 
abundance variations of O, Na, Mg, and Al have been found (Gratton et al. 2001; Ramirez & 
Cohen 2003; Cohen & Melendez 2005). In these unevolved stars, the internal temperatures 
are too low to run the Ne-Na or Mg-Al chains which therefore precludes internal mixing as 
a viable explanation for the star-to-star composition differences. 

In M 4, Smith et al. (2005) found that F varied from star-to-star and that the F abun- 
dance was correlated with O and anticorrelated with Na and Al. Since destruction of F 
is expected to take place during HBB in IM-AGBs (Lattanzio et al. 2004), the observed 
dispersion of F is in qualitative agreement with IM-AGBs being responsible for the glob- 
ular cluster abundance anomalies. However, a quantitative test involving recent yields for 
AGB stars combined with a standard initial mass function showed that the observed abun- 
dance variations cannot be reproduced via pollution from AGB stars (Fenner et al. 2004). 
Denissenkov & Herwig (2003) and Denissenkov & Weiss (2004) also find fiaws in the AGB 
pollution scenario based on calculated yields from AGB models. Ventura & D'Antona (2005) 
caution that theoretical yields from AGB models are critically dependent upon the assumed 
mass-loss rates and treatment of convection such that the predictive power of the current 
AGB models is diminished. That there is still no satisfactory explanation for the star-to- 
star abundance variations seen in every well studied Galactic globular cluster would suggest 
that our understanding of globular cluster chemical evolution and stellar nucleosynthesis is 
incomplete. 

Two neutron-capture elements, rubidium and lead, may offer further clues regarding the 
processes that gave rise to the star-to-star abundance variations and possibly the formation 
of globular clusters. Rb has two stable isotopes, ^^Rb and ^^Rb. While the solar abundance 
of Rb is due to 50% s- process and 50% r- process (Burris et al. 2000), the abundance of 
Rb relative to nearby elements such as Sr, Y, and Zr offers an insight into the neutron 
density at the site of the s-process and therefore the mass of the AGB star due to the 10.7 
yr half-life of ®^Kr (e.g., Tomkin & Lambert 1983, 1999; Lambert et al. 1995; Busso et al. 
1999; Abia et al. 2001). Along the s-process path, Rb is preceded by Kr. The path enters 
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at ^°Kr and exits at either ^^Kr or ^^Kr with ^^Kr providing the branching point. At low 
neutron densities (iV„ < 1 x 10^ cm~^), ^^Kr /9-decays to the stable isotope ®^Rb. At high 
neutron densities, ^^Kr will capture neutrons to form '^^Kr and then *^Kr which /9-decays 
to ^'''Rb (effectively stable with a half-life of 4.7 x 10^° yr). Clearly the isotopic mix of 
Rb depends upon the neutron density. Unfortunately, stellar Rb isotope ratios cannot be 
measured (Lambert & Luck 1976). In the presence of a steady flow along the s-process path, 
the density of a nuclide satisfies the condition aiNi ~ constant, where (jj and Ni are the 
cross-section and abundance of nuclide i respectively. The neutron-capture cross-sections 
differ by a factor of 10 between the two Rb isotopes {asr — (Jgs/lO). The ^^Kr branch does 
not affect the Zr abundances since the low and high neutron density s-process paths converge 
at Sr. Therefore, in a high neutron density environment such as the helium intershell during 
a thermal pulse in IM-AGBs, the Rb abundance may increase by a factor of 10 relative to 
nearby s-process elements such as Sr, Y, and Zr. In reality, the situation is slightly more 
complex since neutron capture on ^"^Kr leads to the ground state of ^^Kr as well as a short 
lived isomeric state that decays to either the ^^Kr ground state or ^^Rb (see Beer & Macklin 
1989 for more details). IM-AGBs of solar metaUicity are expected to have a high neutron 
density with ^^Ne(Q;,n)^^Mg providing the neutron source. (Low-mass AGBs, whose neutron 
source is ^^G(Q;,n)^^0, provide a lower neutron density.) For metal-poor or zero metaUicity 
IM-AGBs, Busso et al. (2001) suggest that such stars do run the s-process though the details 
are model dependent. Nevertheless, the Rb abundance relative to Sr, Y, and Zr (which are 
not affected by the ^^Kr branch) is a potential diagnostic of the s-process site and may offer 
an additional insight into the role of IM-AGBs in the chemical evolution of globular clusters. 

The isotopes of Pb, along with bismuth, comprise the last stable nuclei along the s- 
process path. In low-mass AGB stars, the neutron source is provided by ^"^C whereas in 
IM-AGBs, ^^Ne provides the neutron source with the division occurring at roughly 4 Mq. 
In low-mass AGBs and IM-AGBs of low metaUicity, overabundances of Pb and Bi may be 
expected if the neutron supply per seed exceeds a certain value (e.g., Goriely & Mowlavi 
2000; Travaglio et al. 2001; Busso et al. 2001). Goriely & Siess (2001) suggest that for AGB 
stars with Z < 0.001, the available neutrons per seed nuclei is greater than the number 
required to produce Pb and Bi. Travagho et al. (2001) suggest that metal-poor IM-AGBs 
play only a minor role in the production of Pb though for their 5 Mq model, the Pb yields 
do not change between [Fe/H] = — 1.3 and solar. Herwig (2004) suggest that metal-poor IM- 
AGBs efficiently activate the ^^Ne neutron source though quantitative s-process yields are 
not presented. Busso et al. (2001) predict high yields of Pb from metal-poor IM-AGBs. The 
ratio of Pb/La and Pb/Ba can be used to probe the nature of the s-process in metal-poor 
AGB stars (Gallino et al. 1998; Goriely & Mowlavi 2000). Numerous observational studies 
have found considerable overabundances of Pb in stars that exhibit large s- and/or r-process 
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enhancements (e.g., Cowan et al. 1996; Sneden et al. 2000; Aoki et al. 2000, 2001, 2002; Van 

Eck et al. 2001, 2003; Johnson & Bolte 2002; Lucatello et al 2003; Sivarani et al. 2004; Ivans 
et al. 2005 and references therein). If the globular cluster star-to-star abundance variations 
are due to pollution from metal-poor IM-AGBs, wc may expect large overabundances of Pb 
and a dispersion in Pb abundances despite the absence of variations and excesses in other 
s-process elements. 

In this paper, we present measurements of Rb and Pb in the globular cluster NGC 6752 
as well as measurements of Pb in the globular cluster M 13. While Rb has been measured 
in two globular clusters, uo Gen (Smith ct al. 2000) and NGG 3201 (Gonzalez & Wallcrstein 
1998), as far as we are aware these are the first measurements of Pb in a globular cluster. We 
chose the globular clusters M 13 and NGC 6752 because they exhibit the largest amphtude 
for the Al variation of all the well studied Galactic globular clusters and therefore offer 
the best opportunity to find abundance variations for Rb and Pb. Previous studies of M 13 
include Cohen (1978) and Peterson (1980) who found large Na variations, Shetrone (1996a,b) 
who showed that the Mg isotope ratios were not constant, Cohen & Melendcz (2005) who 
discovered abundance variations in unevolved stars, as well as analyses by Kraft et al. (1992, 
1997), Pilachowski ct al. (1996), and Sneden et al. (2004b). Previous studies of NGG 6752 
include Gottrcll & Da Costa (1981) who discovered the Na and Al enhancements, Suntzeff 
& Smith (1991) who found C and N to systematically vary according to evolutionary status, 
Gratton et al. (2001) and Grundahl et al. (2002) who discovered 0-Al variations in unevolved 
stars, Yong et al. (2003a) (hereafter Y03) who found variations in Mg isotope ratios, Yong 
et al. (2005) (hereafter Y05) who presented evidence for sfight abundance variations of Si, 
Y, Zr, and Ba, and Pasquini et al. (2005) who measured Li in main sequence stars. 



2. Observations and data reduction 

The list of candidates included 5 giants in NGC 6752 previously studied in Y03 and 
Y05, 4 giants in M 13 previously studied by Shetrone (1996a,b), and the comparison star 
HD 141531, a giant whose evolutionary status and stellar parameters are comparable to the 
cluster giants. Though wc were restricted to the brightest giants, the globular cluster stars 
were deliberately selected to span a large range of the star-to-star abundance variations. 
Table 1 contains the list of targets observed using either the Subaru or Magellan telescopes. 

Observations of the M 13 giants and the comparison field star HD 141531 were obtained 
with the Subaru Telescope using the High Dispersion Spectrograph (HDS; Noguchi ct al. 
2002) on 2004 June 1. A 0.4" slit was used providing a resolving power of 90,000 per 4 pixel 
resolution element with wavelength coverage from 4000 A to 6700 A. For the Subaru data. 
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one-dimensional wavelength calibrated normalized spectra were produced in the standard 
way using the IRAF^ package of programs. 

Observations of the NGC 6752 giants were obtained with the Magellan Telescope using 
the Magellan Inamori Kyocera Echelle spectrograph (MIKE; Bernstein et al. 2003) on 2004 
April 3-5. A 0.35" slit was used providing a resolving power of 55,000 in the red and 65,000 

in the blue per 4 pixel resolution element with wavelength coverage from 3800 A to 8500 A. 
While IRAF was used for most of the data reduction, extraction of the Magellan data must 
account for the "tilted" slits, i.e., the lines are tilted with respect to the orders and the tilt 
varies across the CCD. While this is a feature of all cross-dispersed cchcllc spectrographs, 
for MIKE data the tilt is severe. We used the mtools'^ set of tasks written by Jack Baldwin 
to correct for the tilt. Failure to make this correction would result in degradation of the 
spectral resolution as shown in Figure 1. (The magnitude of this effect depends upon the 
aperture size applied to the order being extracted. For the stellar spectra, the decrease in 
spectral resolution would be smaller than for the Th-Ar comparison spectra by roughly a 
factor of 2.) 



3. Analysis 

3.1. Stellair parameters and the iron abundance 

The first step in the analysis was to determine the stellar parameters: the effective tem- 
perature (Teff), the surface gravity (log g), and the microturbulent velocity {^t)- Equivalent 
widths (EWs) were measured for a set of Fcl and Fell lines using routines in IRAF. We 
used the same set of Fe lines presented in Y03. In Figure 2, we compare the measured EWs 
for Fel and Fell lines for the five NGC 6752 giants analyzed in Y03. The EWs measured 
in the Magellan data are in very good agreement with those measured in the VLT data. 
Therefore, for the five NGC 6752 giants, we adopt the same stellar parameters used in Y03. 
In Figure 3, we compare the EWs of Fel and Fell lines for the comparison star HD 141531. 
The EWs measured in the Subaru data are in very good agreement with those measured in 
the VLT data. For HD 141531, we adopt the stellar parameters used in Y03. For the four 
M 13 giants, we determined the stellar parameters using spectroscopic criteria. As in Y03 



^IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under contract 
with the National Science Foundation. 

^http://www.lco.cl/lco/magellan/instruments/MIKE/reductions/mtools.html 
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and Y05, the model atmospheres were taken from the Kurucz (1993) local thermodynamic 
equilibrium (LTE) stellar atmosphere grid and we used the LTE stellar line analysis program 
MoOG (Sneden 1973). For the effective temperature (Teg), we forced the Fel lines to show 
no trend between abundance and lower excitation potential, i.e., excitation equilibrium. To 
set the surface gravity (log(?), wc forced the abundances from Pel and Fell to be equal, i.e., 
ionization equilibrium. We adjusted the microturbulent velocity (^t) until there was no trend 
between abundance and EW. The final [Fe/H] was taken to be the mean of all Fe hnes. Our 
stellar parameters for the M 13 giants compare very well with those derived by Shetrone 
(1996a,b), Sneden et al. (2004b), and Cohen & Melendez (2005). 



3.2. Rubidium abundances 

For the M 13 giants, the HDS spectra did not incorporate the Rb line so we were only 
able to measure Rb in the NGC 6752 giants. The abundances were determined via spectrum 
synthesis of the Rb I line near 7800 A (see Figures 4 and 5). Spectrum synthesis was essential 
for determining accurate abundances due to hyperfine splitting and isotopic shifts as well 
as blending from a stronger Si I line. While the 7800 A Rb I hne is only 3-4% deep relative 
to the continuum, the high quality spectra allow us to measure an abundance from this 
line. Following Tomkin & Lambert (1999), the wavelengths and relative strengths for the 
isotopic and hyperfine structure components were taken from Lambert & Luck (1976) and 
we assumed a solar isotope ratio of ®^Rb/^''Rb = 3. The macroturbulent broadening was 
assumed to have a Gaussian form and was estimated by fitting the profile of the nearby Ni I 
hne at 7798 A. We then generated synthetic spectra and varied the Si and Rb abundances 
to obtain the best match to the observed spectrum. Ideally we would like to measure the Rb 
isotope ratio, ^^Rb/^^Rb, but Lambert & Luck (1976) were unable to measure accurate ratios 
in Arcturus even when using data with superior spectral resolution and signal-to-noise ratio 
(S/N) due to the presence of hyperfine structure and the small isotopic shift. While our tests 
confirmed that we could not measure accurate Rb isotope ratios, we verify the finding by 
Tomkin & Lambert (1999) that the derived Rb abundances arc not sensitive to the assumed 
isotope ratio. Using the Kurucz et al. (1984) solar atlas, we measured an abundance log 
e(Rb) = 2.58 using a model atmosphere with T^fi/logg/it = 5770/4.44/0.85. Our derived 
solar Rb abundance is in very good agreement with the Grevesse & Sauval (1998) value, log 
e(Rb) = 2.60. 

The subordinate Rb i line near 7947 A is weaker by a factor of 2 and is roughly 2% deep 
relative to the continuum. We detect this line in all our spectra and preliminary analyses 
suggest that the abundances derived from this hne agree with those from the 7800 A line (see 
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Figure 6). However, we prefer to restrict our results to the 7800 A line since the 7947 A region 
is more crowded, i.e., the continuum placement strongly affects the derived Rb abundances 
from the 7947 A line. Furthermore, unidentified blends, atmospheric absorption, and fringing 
are more prevalent in this spectral region. (CN lines lie in this region may be absent in these 
metal-poor globular cluster stars.) 

In the subsequent sections, we compare our globular cluster Rb abundances with field 
and cluster stars with [Fe/H] > —2.0 analyzed by other investigators. We now attempt to 
place the various Rb abundance measurements onto a common scale. While the isotope ratio 
of Rb cannot be measured in Arcturus, the elemental abundance ratio is well known. Using 
the Hinkle et al. (2000) Arcturus atlas, we measured an abundance [Rb/H] = —0.55 using 
a model atmosphere with Tf^g /log g/C,t = 4300/1.5/1.55 obtained using the spectroscopic 
criteria described in Section 3.1. The adopted Rb gf value was identical to that used by 
Tomkin & Lambert (1999) and our derived Rb abundance is in very good agreement with 
their measured value, [Rb/H] = —0.58. Since Arcturus is common to both studies and our 
derived abundances are essentially identical, we therefore make no adjustment to the Tomkin 
& Lambert Rb abundances. Gratton & Sncden (1994) use the same gf value, though the 
relative strengths of the hyperfine components differ slightly. We do not adjust their Rb 
abundances. Abia et al. (2001) adopt a Rb gf value that differs from ours, so we adjust 
their Rb abundances by +0.08 dex. It is not clear what Rb gf was used by Gonzalez & 
Wallerstein (1998). Fortunately, Arcturus was also part of their sample. They derived an 
abundance [Rb/H] = —0.45 and so we adjust their Rb abundances by —0.1 dex. Smith et al. 
(2000) find [Rb/H] = —0.52 and so we do not adjust their Rb abundances. 



3.3. Lead abundances 

The Pb abundances were determined via spectrum synthesis of the Pbl line near 
4058 A (see Figures 7 and 8). Abundances from the Pb hne near 3683 A could not be 
determined due to the lack of flux in the blue for these cool giants. While the region cen- 
tered near 4058 A is crowded with molecular hues of CH as well as atomic hues from Mg, 

Ti, Mn, Fc, and Co, our syntheses provide a very good fit to the region demonstrating that 
reliable Pb abundances can be extracted. The macroturbulent broadening was estimated by 
fitting the profiles of the nearby lines. We adopted the same gf value used by Aoki et al. 
(2000, 2001, 2002). Following Aoki et al. (2002), our synthesis accounted for the hyperfine 
and isotopic splitting as well as the isotopic shifts. (The stable isotopes are ^°''Pb, ^°^Pb, 
^°^Pb, and ^'^^Pb.) We again assumed a solar isotope ratio for Pb though as with Rb, our 
tests confirmed that the derived elemental Pb abundance was not sensitive to this choice. 
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For the solar Pb abundance, we adopted log e(Pb) = 1.95 from Grevesse & Sauval (1998). 

In the subsequent sections, we compare the globular cluster Pb abundances with those 
obtained in field stars with [Fe/H] > —2.0 by other investigators. However, we further 
restrict the comparison by avoiding stars with known s-process enhancements leaving only 
a handful of stars from two studies, Sneden et al. (1998) and Travaglio et al. (2001). Again 

we attempt to put the Pb abundance measurements onto a common scale. Our gf value 
is identical to that used by Sncdcn ct al. (1998) so we do not adjust their Pb abundances. 
Travaglio et al. (2001) used a different Pb line (3683 A) and they did not list the adopted 
gf value. Since there arc no stars common to both analyses, we do not make an adjustment 
to their Pb abundances and caution that there may be a systematic offset. 

3.4. Additional elements 

We also measured abundances for Al, Si, Y, Zr, La, and Eu in the M 13 giants and the 
comparison star HD 141531 using the same lines presented in Y05. These measurements 
were performed to ensure that the abundances would be on the same scale as the NGC 6752 
giants studied in Y05. Zr was chosen because we compare the Rb and Zr abundances to 
look for a large ratio [Rb/Zr] as well as a detectable dispersion, i.e., the hallmark of a high 
neutron-density environment and the possible signature of pollution from IM-AGBs. Al, Si, 
and Y were also chosen because Y05 found evidence for correlations between Al and Si, Al 
and Y, and Al and Zr in NGC 6752. While our sample size in M 13 is small, it would be 
interesting to see if similar correlations are present. La and Eu were measured since these 
neutron-capture elements offer an insight into the ratio of s-process to r-process material. 
Furthermore, the ratio [Pb/La] has been used to test predictions from AGB models. In 
Table 2 we present our measured elemental abundances for Al, Si, Rb, Y, Zr, La, Eu, and 
Pb in the program stars. The adopted solar abundances for Al, Si, Y, Zr, La, and Eu were 
6.47, 7.55, 2.24, 2.60, 1.13, and 0.52 respectively 

We attempt once more to put the abundance measurements onto a common scale by 
considering the gf values used by the various studies to which we compare our abundances. 
The element we focus upon is Zr (in order to compare [Rb/Zr] between the samples). We 
shift all the Zr abundances onto the Smith et al. (2000) scale in order to compare with their 
theoretical predictions for [Rb/Zr] from low and intermediate-mass AGBs (their Figure 14). 
(Our Rb abundances were already on the Smith scale.) Abundance measurements for Zr 
are complicated by the fact that the laboratory Zr gf values from Biemont et al. (1981) are 
smaller than the solar (?/- values by 0.41 dex (Tomkin & Lambert 1999). We must therefore 
take care and account for both the adopted solar abundance and the gf value. Smith et al. 
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(2000) adopt the Biemont et al. (1981) gi/-values and a solar abundance log e0(Zr) = 2.90. 
We used the Biemont et al. (1981) g^/- values and a solar abundance log e0(Zr) = 2.60 
(Grevesse & Sauval 1998). Therefore we adjust our Zr abundances by —0.30 dex to ensure 
that we are on the same scale as Smith et al. (2000). Similarly, we adjust the Zr abundances 
of Gratton & Sneden (1994) and Abia et al. (2001) by —0.30 dex since they adopt the same 
gf values and a very similar solar abundance used in our analysis. Gonzalez & Wallerstein 
(1998) adopt the Grevesse & Sauval (1998) solar abundance and a different gf value so we 
adjust their Zr abundances by +0.06 dex. Tomkin & Lambert (1999) adopt the Grevesse 
& Sauval (1998) solar abundance and a different gf value so we adjust their Zr abundances 
by +0.11 dex. These adjustments are substantial. When we compare the ratio [Rb/Zr], we 
also consider how the comparison would fare had we not made these abundance corrections. 
To assess the validity of these adjustments, we measured the Zr abundance for Arcturus and 
found [Zr/H] = —0.66. Smith et al. measured [Zr/H] = —0.96, Gonzalez & Wallerstein found 
— 1.18, and Tomkin & Lambert found —1.00. Therefore, applying the abundance corrections 
based on the gf values and solar abundances ensures that the Zr abundances are on the 
Smith et al. (2000) scale, e.g., for Arcturus we find —0.96 (this study), —1.12 (Gonzalez & 
Wallerstein), —0.89 (Tomkin & Lambert), and —0.96 (Smith et al.). Interestingly, if we had 
used the same solar abundance as Smith et al. (2000), our [Zr/Fe] abundances in Y05 would 
have been closer to the [Y/Fe] values and for M 13 we would have found [Y/Fe] ~ [Zr/Fe]. 
In Figure 9, we plot [Zr/Fe] versus [Fe/H] with the abundances shifted to the Smith et al. 
(2000) scale. At the metallicity of NGC 6752 and M 13, the field and cluster stars have 
similar ratios of [Zr/Fe]. The comparison field star HD 141531 has [Zr/Fe] almost identical 
to the globular clusters. 

As in Y03, we estimate the internal errors in the stellar parameters to be Tcs ± 50K, 
logg ± 0.2, and =t 0.2. In Table 3, we show the abundance dependences upon the model 
parameters. 

4. Discussion 
4.1. Rubidium 

While our mean Rb abundance for NGC 6752 is [Rb/Fe] = -0.17 ± 0.06 (a = 0.14), 
the abundances appear to concentrate around two distinct values. There are two stars with 
[Rb/Fe] ~ -0.02 and three stars with [Rb/Fe] ~ -0.25. The two stars with the higher Rb 
abundances do not have the highest Al abundances and the three stars with the lower Rb 
abundances are not exclusively the stars with the lowest Al abundances. Given the weakness 
of the Rb line, the uncertainties in the derived Rb abundances (see Table 3), and the small 
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sample size, it is unlikely that the Rb abundances show a dispersion in NGC 6752. Nor do 
we find evidence for a correlation between [Al/Fe] and [Rb/Fc], though we recognize that our 
sample size (5 stars) is much more limited than in Y03 and Y05 (38 stars). Unfortunately, 
observations of M 13 and the comparison field star HD 141531 did not incorporate the Rb 
line. 

In Figure 10, we compare our Rb abundances with those measured in dwarfs and giants 
in the disk and halo (Gratton & Sneden 1994 and Tomkin & Lambert 1999), globular cluster 
giants (Gonzalez & Wallerstein 1998 and Smith et al. 2000), and carbon-rich AGB stars 
(Abia ct al. 2001). (While wc retain their stars in the plots, wc note that the Abia sample 
contains very different objects that are difficult to analyze compared to the dwarfs and 
giants considered in the other studies.) Recall that we have made small adjustments to the 
Rb abundances in an attempt to place them onto a common scale. At the metallicity of 
NGC 6752 ([Fe/H] = -1.6), our two stars in NGC 6752 with the highest [Rb/Fe] ratios have 
abundances compatible with the lower envelope of the Tomkin & Lambert (1999) sample. 
The two NGC 6752 stars also exhibit very similar abundances [Rb/Fc] to the Gonzalez & 
Wallerstein (1998) and Smith et al. (2000) globular cluster giants. Our three stars with the 
lower [Rb/Fe] ratios appear unusual compared to the Tomkin & Lambert sample. Only 1 
star in the Abia et al. (2001) sample has [Fe/H] < —1.0 and it is interesting that it has an 
abundance ratio [Rb/Fe] similar to those measured in NGC 6752. In general, Rb is not a 
well studied element and the comparison data are limited. 

For NGC 6752, we find a mean abundance [Rb/Zr] = -0.12 ± 0.06 {a = 0.13). (This 
abundance ratio has been shifted to the Smith et al. (2000) scale.) For the five NGC 6752 
giants, the ratio [Rb/Zr] appears to show a dispersion. We suspect that this is attributable 
to measurement uncertainties (primarily for Rb) rather than reflecting a real star-to-star 
scatter. We do not find a correlation between [Al/Fe] and [Rb/Zr]. In Figure 11, we compare 
the abundance ratio [Rb/Zr] between NGC 6752 and various field and cluster stars. Note 
that in this Figure we have shifted all abundances onto the Smith et al. (2000) scale since 
we will utilize their theoretical predictions from low and intermediate-mass AGBs. At the 
metallicity of NGC 6752, we find that the two stars in NGC 6752 with the highest values 
of [Rb/Fe] also have the highest values of [Rb/Zr]. These two stars have similar [Rb/Zr] 
ratios to the Gratton & Sneden (1994) and Tomkin & Lambert (1999) samples at the same 
metaUicity. While the u; Cen giants have abundance ratios [Rb/Zr] shghtly lower than NGC 
6752, this time the NGC 3201 giants appear to have much higher ratios of [Rb/Zr]. Note 
that the [Rb/Zr] ratios in NGC 3201 appear similar to the highest values seen in the Tomkin 
& Lambert (1999) sample. Unfortunately, the only star in the Abia et al. (2001) sample 
with [Fe/H] < —1.0 does not have a Zr measurement. However, it does have [Rb/Sr] = —0.5 
and [Rb/Y] = -0.6. If we assume for this star [Rb/Zr] = < [Rb/Sr], [Rb/Y]> = -0.55, then 
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the abundance is much lower than NGC 6752. 

For elements heavier than Si, globular clusters and field stars tend to have very similar 
abundance ratios [X/Fe] at a given [Fe/H] (Gratton et al. 2004; Sneden et al. 2004a). Al- 
though the scatter is large and the sample sizes are limited, it would appear that cluster stars 
probably have similar, or perhaps slightly lower abundance ratios of [Rb/Fe] and [Rb/Zr] 
compared to field stars at a given [Fe/H]. 

Recall that we made substantial adjustments to the Zr abundance. While consideration 
of the Arcturus Zr abundances would appear to validate this adjustment, we briefiy consider 
how the comparison of [Rb/Zr] would have fared if these corrections were not applied. In 
this case, the ratio [Rb/Zr] would decrease by roughly 0.3 dex for NGC 6752 as well as 
for the Gratton & Sneden (1994) and Abia et al. (2001) samples. The Tomkin & Lambert 
(1999) and Gonzalez & Wallerstein (1998) samples would increase by roughly 0.15 dex. NGC 
6752 would therefore have unusually low ratios [Rb/Zr] compared to field stars at the same 
metallicity. Similarly, the cu Cen compositions would be unusually low though comparable 
to NGC 6752. The NGC 3201 giants would then have very high [Rb/Zr] ratios compared to 
other globular clusters and field stars at the same metaUicity. NGC 3201 is pecuhar since 
it has a retrograde orbit and may have been a captured cluster (van den Bergh 1993). The 
capture hypothesis could not be demonstrated by Gonzalez & Wallerstein (1998) who found 
no unusual abundance ratios. 

Smith et al. (2000) compare [Rb/Zr] in u; Cen with predictions from AGB models 
with various initial masses and initial metallicities (their Figure 14). Their Figure clearly 
shows how the ratio [Rb/Zr] can vary by nearly a factor of 10 depending on whether a 
low-mass (1.5Mo) or high-mass (5M0) AGB model is synthesizing the s-process elements. 
As anticipated from the arguments given in Section 1, high- mass AGB models produce high 
[Rb/Zr] while low-mass AGB models produce low [Rb/Zr]. The magnitude of the difference 
in the predicted [Rb/Zr] between low- and high-mass AGB models does not significantly 
change as the metallicity decreases from [Fe/H] = —0.5 to [Fe/H] = —2.0. Comparing the 
observed abundances with the model predictions in Smith et al. reveals that low-mass AGB 
stars (1-3 M0) are responsible for the synthesis of the s-process elements in u; Cen. Inspection 
of Figure 14 in Smith et al. (2000) shows that at the metaUicity of NGC 6752, [Fe/H] = 
—1.6, our measured ratio [Rb/Zr] = —0.12 is compatible with the s-process elements being 
synthesized in low-mass AGB stars though the assumed mass of the ^'^C pocket is critical. 
Predictions assuming a standard treatment for the ^^C pocket or the ^^C pocket increased 
by a factor of 2 both suggest AGB stars with < 3 Mq are responsible for the [Rb/Zr] ratios 
seen in NGC 6752. When the ^^C pocket is diminished by a factor of 3, the AGB stars with 
masses > 3 Mq may explain the observed [Rb/Zr]. When we return our Zr abundances to 
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the original scale, [Rb/Zr] = —0.42, the ratio in NGC 6752 is only compatible with low-mass 
AGB stars. We note that the highest values of [Rb/Zr] seen in the Gratton & Sneden (1994), 
Gonzalez & Wallerstein (1998), and Tomkin & Lambert (1999) samples all greatly exceed 
the 5 M0 AGB model predictions. Such a discrepancy serves as a useful reminder of the 
unfortunate reality that the detailed yields of s-process elements from AGB stars may be 
very model dependent (Busso et al. 2001; Ventura & D'Antona 2005). 

4.2. Lead 

In NGC 6752, the mean Pb abundance is [Pb/Fe] = -0.17 ± 0.04 {a = 0.08) and 
in M 13 the mean abundance is [Pb/Fe] = —0.28 ± 0.03 {a = 0.06). Given the fairly 
large measurement uncertainty for Pb (sec Tabic 3), neither NGC 6752 nor M 13 show any 
evidence for a dispersion in Pb abundances, though our sample sizes for both clusters are 
small. Furthermore, the [Pb/Fe] ratios are very similar for these two clusters. As with Rb, 
there is no evidence that the ratio [Pb/Fe] is correlated with [Al/Fe]. We note that one star, 
NGC 6752 PDl, has lower ratios of both [Rb/Fe] and [Pb/Fe] relative to other giants in this 
cluster. This subtle composition difference probably arises from uncertainties in the stellar 
parameters rather than representing a genuine difference. The comparison field giant HD 
141531 has a ratio [Pb/Fe] essentially identical to the globular cluster giants. 

In Figure 12, we compare our Pb abundances with values measured by Sneden et al. 
(1998) and Travagho et al. (2001). While Pb has been measured in numerous stars with 
large s-process enhancements, it has been largely neglected in normal field stars presumably 
due to the difficulty of the measurement. For HD 126238, the Pb abundance measured by 
Sneden et al. (1998) is very similar to the globular cluster giants and HD 141531. The Pb 
abundances measured by Travaglio et al. (2001) in field stars are larger than those measured 
in the globular clusters. The star with [Pb/Fe] = 0.6 is a CH star with excess C and Ba and 
should not be considered a normal field star. Aside from the CH star, there are three stars 
with upper limits and another three Pb detections. Recall that there are no stars common 
to both studies and that the Pb gf value was not published. Travaglio et al. (2001) suggest 
that some of the Pb detections may be uncertain and therefore, the offset between the Pb 
abundances may be due to measurement errors and/or the adopted gf value. 

Travaglio et al. (2001) not only measured Pb abundances in a handful of stars, but they 
calculated the Galactic chemical evolution of Pb from a detailed model. In their Figure 4, 
they plot the expected run of [Pb/Fe] versus [Fc/H] for the halo, thick disk, and thin disk. 
Since their prediction integrates over all AGB masses, it would be useful to learn how the 
predicted curve would differ (if at all) if the calculation was performed using low-mass or 
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high-mass AGB models exclusively as Smith et al. (2000) have done. At [Fe/H] = —1.6, the 
Travaglio et al. (2001) model predicts an abundance ratio [Pb/Fc] ~ —0.1. This prediction 
is in very good agreement with the values measured in M 13, NGC 6752, HD 126238, and 
HD 141531. This agreement may be regarded as evidence that globular cluster stars have 
virtually identical Pb abundances as normal field stars. 

In normal field stars, Pb has been less studied than Rb. Clearly, it would be of great 
interest to have additional Pb measurements in field and cluster stars. In Figures 7 and 8, 
our syntheses indicate that for cool giants in the metallicity regime —2.0 < [Fc/H] < —1.0, 
reliable Pb abundances can be measured even in stars that do not have large Pb or s-process 
enhancements. 

The ratio [Pb/La] may offer further clues regarding the nature of the s-process in the 
AGB stars. Van Eck et al. (2003) found some stars with ratios of [Pb/La] > +1.5, in 
agreement with predictions from metal-poor AGB models (Gallino et al. 1998; Goriely & 
Mowlavi 2000). However, Aoki et al. (2002) and Van Eck et al. (2003) also found a large 
spread in the ratio [Pb/Ba]. In some stars, the ratio [Pb/Ba] was sub-solar. Our mean ratio 
[Pb/La] for M 13 is -0.36 ± 0.05 (cr = 0.10). For NGC 6572, our mean ratio [Pb/La] = 
—0.23 ± 0.04 ((7 = 0.09). In both clusters, the mean ratios are similar and we note that they 
are both sub-solar and comparable to the lowest ratios seen in the Van Eck et al. (2003) 
sample. Curiously the subsample in Van Eck et al. (2003) with [Pb/La] < had extreme 
enhancements for [Pb/Fe] and [La/Fc]. The comparison field star HD 141531, has [Pb/La] 
= —0.20 which is similar to the value seen in the globular clusters. 

4.3. Additional elements 

While our sample in M 13 consists of only 4 stars, they span the extremities of the Al 
variation. As in Y05, we again find that the most Al-rich stars may also exhibit slightly 
higher Si abundances than the most Al-poor stars. Further measurements of Al and Si in 
a large sample of stars in M 13 would be of great interest to verify whether the correlation 
between Al and Si seen in NGC 6752 (Y05) is also present in M 13. The correlations between 
Al and Y as well as Al and Zr found in NGC 6752 do not appear to be present in the small 
M 13 sample. Cohen & Melendez (2005) measured abundances in 25 stars in M 13 from the 
main sequence turn-off to the tip of the RGB. They were unable to measure Al in most stars. 
When we consider their derived abundances, there appears to be an anticorrelation between 
O and Si as well as O and Y. Though the anticorrelation is driven primarily by the one star 
with unusually low [0/Fe], such trends are intriguing and warrant further investigation. 
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As noted in previous investigations of these clusters, the ratio of s-process to r-process 
material, [La/Eu], is sub-solar but greater than the scaled solar pure r-process value. For 
NGC 6752 and M 13, the observed ratios of [La/Eu] show that AGBs have contributed to 
their chemical evolution. The ratio of La/Eu in HD 141531 again confirms that it is a normal 
field star. 



4.4. Consequences for the IM-AGB pollution scenario 

In Y03, we measured Mg isotope ratios in bright giants in NGC 6752. We found that 
the ratio varied from star-to-star. Specifically, ^'^Mg was anticorrelated with Al, ^^Mg was 
correlated with Al, and ^^Mg was not correlated with Al. As previously seen by Shetrone 
(1996b) in M 13, these isotope ratios reveal that the Al enhancements result from proton 
capture on the abundant ^"^Mg. Proton capture on ^"^Mg within the Mg-Al chain is predicted 
to only occur in AGB stars of the highest mass at their maximum luminosity (Karakas & 
Lattanzio 2003) . So we suggested that the abundance variations were due to differing degrees 
of pollution by IM-AGBs, an idea originally proposed by Cottrell & Da Costa (1981). These 
IM-AGBs must have the same iron abundance as the present generation of cluster stars 
otherwise there would also be a star-to-star abundance variation of Fe. (The same argument 
applies to whatever stars are believed to be the source of the pollutants.) 

The Mg isotope ratios presented in Y03 offered further clues to globular cluster chemical 
evolution. At one extreme of the abundance variation are cluster stars with O, Na, Mg, and 
Al compositions in accord with field stars at the same metallicity. We called such stars 
"normal" in anticipation that proton capture nucleosynthesis can produce 0-poor, Na-rich, 
Mg-poor, and Al-rich material. At the other extreme of the abundance variations are the 
stars with high Na, high Al, low O, and low Mg. We referred to these stars as "polluted". 
The pollution may have occurred via either the evolutionary or primordial scenario. In 
"normal" stars, we found ratios ^^Mg/^^Mg and ^^Mg/^^Mg that exceeded field stars at the 
same metaUicity (Yong et al. 2003b). Of equal importance was the fact that these isotope 
ratios greatly exceeded predictions from metal-poor supernovae. We therefore suggested that 
these unusually high isotope ratios could be explained if a previous generation of IM-AGBs 
of the highest mass polluted the natal cloud from which the cluster formed. The ejecta from 
this previous generation must have been thoroughly mixed before the present generation 
of stars began to form. This previous generation of IM-AGBs are probably responsible for 
much of the Na, Al, and N as well as ^^Mg and ^^Mg. 

Our working hypothesis is that IM-AGBs played two crucial roles in globular cluster 
chemical evolution. Firstly, a prior generation of very metal-poor IM-AGBs are required to 
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produce the high ^^Mg/^^Mg and ^^Mg/^^Mg seen in "normal" stars. Secondly, a generation 
of IM-AGBs with the same Fe abundance as the present cluster members pollutes the cluster 
environment. Differing degrees of pollution of natal clouds then produce the star-to-star 
abundance variations. The dispersion in the F abundances and the correlation between F 
and O in M 4 (Smith et al. 2005) appear to confirm the role of IM-AGBs in producing 
the abundance variations. However, not all the abundance patterns observed in globular 
clusters can be matched by the current theoretical yields from IM-AGBs (Denissenkov & 
Herwig 2003; Denissenkov & Weiss 2004) nor can the abundance patterns be reproduced by 
chemical evolution models (Fenner et al. 2004). 

From a qualitative viewpoint, metal-poor IM-AGBs may produce s-process elements 
via the ^^Ne neutron source. If activated, the ^^Ne neutron source produces large amounts 
of Rb/Zr due to a critical branching point at ^^Kr as described earher. Theoretical models 
by Busso et al. (2001) suggest that metal-poor IM-AGBs do run the s-process though the 
specific yields depend on the details. 

Similarly, a quafitative assessment suggests that metal-poor IM-AGBs will produce Bi 
and Pb if the neutrons per seed nuclei exceed a certain value. In this case, Bi and Pb 
may show large enhancements with other s-process elements showing only modest overabun- 
dances. Again, theoretical models can be found in which metal-poor IM-AGBs do produce 
lead (e.g., Goriely & Siess 2001 and Busso et al. 2001). 

In M 13 and NGC 6752, we did not find high ratios of [Rb/Fe], [Rb/Zr], or [Pb/Fe] 
compared with field stars at the same metallicity. If metal-poor IM-AGBs are responsible 

for the globular cluster star-to-star abundance variations, then our measurements strongly 
suggest that such stars do not synthesize significant quantities of Rb or Pb. Also, if metal- 
poor IM-AGBs are responsible for the large abundances of ^^Mg and ^^Mg in "normal" 
cluster stars, then they do not synthesize Rb or Pb. Alternatively, if metal-poor IM-AGBs 
do synthesize significant quantities of Rb and Pb, then they cannot be responsible for the 
abundance anomalies seen in globular clusters. 

Of course the possibility remains that the predicted yields of Rb and Pb from IM- 
AGBs are unreliable and/or model dependent (e.g., Ventura & D'Antona 2005). It has been 
suggested that metal-poor IM-AGBs will not produce any s-process elements. As the mass of 
the AGB star increases, the size of the He intershell region decreases as does the duration of 
the thermal pulse (Lattanzio et al. 2004). Detailed predictions of the yields from metal-poor 
IM-AGBs by independent groups are required. Our observed Rb and Pb abundances may 
serve to constrain these models. 

Finally, we note that Rb and Zr are synthesized in IM-AGBs as well as via the weak 
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S' process, i.e., He core burning in massive stars. However, Travaglio et al. (2004) suggest 
that the weak s-process does not contribute to Zr but a hghter element primary process in 
massive stars may be responsible for up to 18% of the solar abundance of Zr. Chieffi & 
Limongi (2004) present detailed yields from massive stars for a range of metallicities and 
masses. At Z = 0, massive stars are predicted to produce low ratios, [Rb/Zr] = —0.5 to —2.4. 
If IM-AGBs have contributed to the chemical evolution of NGC 6752, they would increase 
the ratios of [Rb/Zr] above those produced by the supernovae. If IM-AGBs have not played a 
role in the chemical evolution, perhaps we can use [Rb/Zr] to probe the mass and metallicity 
range of the previous generation of massive stars. For Z — 10"^, the predicted yields are 
independent of mass with [Rb/Zr] = —0.23. However, we note that for other metallicities, 
there is a mass- metallicity degeneracy for the [Rb/Zr] yields that limits their use in probing 
the previous generation of supernovae. 



5. Concluding remarks 

We show for the first time the uniformity of the neutron-capture elements Rb and Pb in 
NGC 6752 and M13, the two globular clusters that exhibit the largest dispersion for Al. We 
also find the ratio [Rb/Zr] to be constant. None of the abundance ratios [Rb/Fe], [Rb/Zr], 
or [Pb/Fe] are correlated with [Al/Fe] and the Rb and Pb abundances show sub-solar ratios 
[X/Fe]. If metal-poor IM-AGBs produce large amounts of Pb and Rb as well as high ratios 
of [Rb/Zr], then such stars are not responsible for the abundance variations, a conclusion 
already suggested by Denissenkov & Herwig (2003), Denissenkov & Weiss (2004), and Fenner 
et al. (2004). If metal-poor IM-AGBs are responsible for the abundance variations, then they 
cannot produce overabundances of Rb or Pb. 

For elements heavier than Al, previous studies have shown that field and cluster stars 
generally have the same abundance ratios [X/Fe] at a given [Fe/H]. While our sample size is 
small and the data for comparison field stars are limited, the two clusters we have studied 
have Rb abundance ratios [Rb/Fe] and [Rb/Zr] in reasonable agreement with the general 
field population (the clusters may have shghtly lower ratios). The Pb abundance ratio 
[Pb/Fe] in globular clusters is in very good agreement with the limited sample of field stars. 
At the metallicity of M 13 and NGC 6752, their Pb abundances are well matched by the 
predictions from the chemical evolution model by Travaglio et al. (2001). In order to further 
our understanding of stellar nucleosynthesis and the chemical evolution of field and cluster 
stars, additional measurements of Rb and Pb in normal stars and globular clusters are 
welcomed as are further theoretical efforts to calculate the yields from metal-poor IM-AGBs. 
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Fig. 1. — Profiles of extracted TliAr lines from the same exposure taken with Magellan-MIKE 
scaled to the same peak intensity. The solid black line shows the ThAr lines extracted without 
taking into account the tilted slits. The dashed red line shows the ThAr lines extracted using 
the MTOOLS package. Note how the FWHM decreases when the correction is performed. 
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Fig. 2. — Comparison of the equivalent widths of Fe I and Fe II lines between this study 
(Magellan data) and Y03 (VLT data) for the sample of NGC 6752 giants. 
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Fig. 3. — Comparison of the equivalent widths of Fe I and Fe II hnes between this study 
(Subaru data) and Y03 (VLT data) for HD 141531. 
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Fig. 4. — Spectrum of NGC 6752 B702 near the 7800 A Rb feature. Synthetic spectra with 
differing Rb abundances are shown. 
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Fig. 5.— Same as Figure 4 but for star NGC 6752 B1630. 
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Fig. 6. — Spectrum of NGC 6752 B708 near the 7947 A Rb feature. Synthetic spectra with 
differing Rb abundances are shown. The lower panel is identical to the upper panel except 
for the y-range. 
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Fig. 7. — Spectrum of NGC 6752 B708 near the 4058 A Pb feature. Synthetic spectra with 
differing Pb abundances are shown. 



- 30 - 




0.2 - 



4056 



4057 4058 
Wavelength(A) 



4059 



Fig. 8. — Same as Figure 7 but for star M 13 L973. 
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Fig. 9. — [Zr/Fe] versus [Fe/H]. Closed black circles show our measurements for NGC 6752, 
the open black circles show M 13, the open black square shows HD 141531, the closed green 
triangles are from Gratton Sz Sneden (1994), the green asterisks represent data from Gonzalez 
& Wallerstein (1998), open blue triangles are from Tomkin & Lambert (1999), red plus signs 
represent data from Smith et al. (2000), and filled red squares show data from Abia et al. 
(2001). A representative error bar is shown and the Zr abundances have been shifted onto 
the Smith et al. scale (see text for details). 
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Fig. 10. — [Rb/Fe] versus [Fe/H]. Closed black circles show our measurements for NGC 6752, 
the closed green triangles are from Gratton & Sneden (1994), the green asterisks represent 
data from Gonzalez & Wallerstein (1998), open blue triangles are from Tomkin & Lambert 
(1999), red plus signs represent data from Smith et al. (2000), and filled red squares show 
data from Abia et al. (2001). A representative error bar is shown. 
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Fig. 11. — [Rb/Zr] versus [Fe/H]. The symbols are the same as in Figure 10. Note that we 
shifted all abundances onto the Smith et al. (2000) scale (see text for details). 
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Fig. 12. — [Pb/Fe] versus [Fe/H]. Closed black circles, open blue triangles, and the closed 
red square show our measurements for NGC 6752, M 13, and the comparison field star HD 
141531. The closed green triangles and green arrows (upper limits) are from Travaglio et al. 
(2001) while the purple asterisk shows HD 126238 from Sneden et al. (1998). A representative 
error bar is shown. 
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Table 1. Exposure times and stellar parameters. 



Star 


Telescope 


Exposure 


S/N*^ 


S/N*^ 




lo 


g 9 


6 


[Fe/H] 






Time (mill) 


4()5()A 


78()()A 


K 






km 




M 13 L598 


Subaru 


60 


39 




3900 


0. 


,00 


2.25 


-1.56 


M 13 L629 


Subaru 


50 


37 




3950 


0. 


,20 


2.25 


-1.63 


M 13 LTO*' 


Subaru 


60 


39 




3950 


0, 


,30 


2.25 


-1.59 


M 13 L973'= 


Subaru 


64 


36 




3920 


0. 


,30 


2.35 


-1.61 


NGC 6752 B702 


Magellan 


20 


42 


197 


4050 


0. 


,50 


2.10 


-1.58 


NGC 6752 B708 


Magellan 


50 


55 


402 


4050 


0. 


,25 


2.20 


-1.63 


NGC 6752 PDl 


Magellan 


26 


47 


197 


3928 


0, 


,26 


2.70 


-1.62 


NGC 6752 B1630 


Magellan 


27 


44 


247 


3900 


0. 


,24 


2.70 


-1.60 


NGC 6752 B3589 


Magellan 


21 


46 


240 


3894 


0. 


,33 


2.50 


-1.59 


HD 141531 


Subaru 


10 


85 




4273 


0. 


,80 


1.90 


-1.72 



^S/N values arc per pixel. 
''Alternative name 11-67. 
■^Alternative name 1-48. 



Table 2. Elemental abundances. 



Star 


[Al/Fe] 


[Si/Fe] 


[Rb/Fe] 


[Y/Fe] 


[Zv/Fe]" 


[Zr/Fe]'' 


[La/Fe] 


[Eu/Fe] 


[Pb/Fe] 


M 13 L598 


0.24 


0.23 




-0.10 


0.27 


-0.03 


0.00 


0.40 


-0.30 


M 13 L629 


0.74 


0.31 




-0.05 


0.25 


-0.05 


0.06 


0.44 


-0.20 


M 13 L70 


1.27 


0.33 




-0.12 


0.24 


-0.06 


0.14 


0.43 


-0.27 


M 13 L973 


1.28 


0.35 




-0.05 


0.17 


-0.13 


0.12 


0.52 


-0.35 


NGC 6752 B702 


1.04 


0.43 


-0.02 


0.04 


0.32 


0.02 


0.09 


0.28 


-0.10 


NGC 6752 B708 


1.23 


0.35 


-0.22 


0.04 


0.30 


0.00 


0.03 


0.28 


-0.15 


NGC 6752 PDl 


1.08 


0.38 


-0.25 


0.07 


0.23 


-0.07 


0.07 


0.32 


-0.30 


NGC 6752 B1630 


0.82 


0.41 


-0.03 


0.04 


0.19 


-0.11 


0.05 


0.34 


-0.10 


NGC 6752 B3589 


0.77 


0.43 


-0.32 


0.10 


0.20 


-0.10 


0.08 


0.35 


-0.18 


HD 141531 


0.01 


0.23 




-0.13 


0.14 


-0.16 


0.00 


0.33 


-0.20 



^Zr abundances using the gf values and solar abundance assumed in Y05. 
''Zr abundances when shifted onto the Smith et al. (2000) scale. 
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Table 3. Abundance dependences on model parameters for NGC 6752 B1630. 



Species 


Teff + 50K 


log 5 + 0.2cgs 


£,t + 0.2km 


[Fe/H] 


0.02 


0.01 


-0.03 


[Al/Fe] 


0.05 


-0.01 


-0.01 


[Si/Fe] 


-0.01 


0.03 


-0.01 


[Rb/Fe] 


0.09 


0.03 


-0.01 


[Y/Fe] 


0.01 


0.06 


-0.04 


[Zr/Fe] 


0.13 


0.02 


-0.01 


[La/Fe] 


-0.02 


0.05 


-0.05 


[Eu/Fe] 


-0.04 


0.05 


-0.04 


[Pb/Fe] 


0.12 


-0.05 


-0.02 



